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The role of the orbital degrees of freedom is studied theoret- 
ically for the spin dynamics of Ri- x A x Mn03. Based on the 
meanfield solution, an RPA calculation has been done and 
it is found that the d x2 _ y2 -type orbital is essential for the 
double-exchange (DE) interactions, i.e., the DE is basically 
two-dimensional interaction. Based on this results compared 
with experiments, we propose that the orbital wavefunction is 
d a; 2_ 1( 2-type locally even in the metallic ferromagnetic state, 
which fluctuate quantum mechanically. Well agreement of 
the estimation with experiments suggest that the Jahn- Teller 
phonon has less importance on the spin dynamics. 
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Doped manganites Ri- x A x Mn03 (i?=La, Pr, Nd, Sm 
; A= Ca, Sr, Ba) have recently attracted considerable 
interests due to the colossal magnetoresistance (CMR) 
observed near the ferromagnetic (spin _F-type) transi- 
tion temperature T c Er1j. It is now recognized that the 
most fundamental interaction in these materials is the 
double exchange-interaction (DE), which connects the 
transport and magnetismd. Therefore the magnetism 
is a key issue to reveal the mechanism of CMR. Es- 
pecially, rich magnetic phase diagrams have been clar- 
ified over the wide range of the concentration x and also 
the bandwidth. With increasing x, the parent insula- 
tor with a layered antiferromagnetism (spin A-type AF) 
changes into a ferromagnetic metal (FM).u In addition to 
these well-known, magnetic phases,.-spin A-type AF lin 
La 1 _- e Sr- e Mn0 3 ,El Pri_ x Sr x Mn0 3 ,§B Ndi_ x Sr a Mn0 3 &9) 
and the rod type antiferromagnetism (spin C-type AF, 
in Ndi-xSr^MnOalilj) were recently found in the moder- 
ately doped metallic region (0.5 < x < 0.8). 

The neutron scattering experiments have revealed the 
spin wave excitation at low temperatures, which de- 
pends sensitively on the doping x and the magnetic 
structuraijil3. In the spin A-type AF for small x, the dis- 
persion is two dimensional while it becomes isotropic in 
FM state. The spin stiffness, however, stays almost con- 
stant up to a: = 0.125 where the phase transition between, 
the insulating and metallic ferromagnetic states occurs .Ell 
This phase transition is accompanied with that of the or- 
bital structureO In the FM state, the orbital ordering 
disappears and the spin stiffness begins to increase. In 
the metallic A- type AF (AFM) state for .higher x, the 
dispersion becomes again two dimensional £-3 In this pa- 
per we report the calculation of the spin wave dispersion 
by changing x and taking into account the orbital struc- 
ture. The calculated x-dependence of the spin stiffness 



agrees quantitatively with the observed one for x > 0.2 
where the double-exchange interaction dominates. This 
x-dependence strongly supports the large orbital polar- 
ization, which is d x 2_ y 2 at least locally. Therefore the 
double-exchange interaction is basically two-dimensional. 

We previously reported a meanfield theory (MFT) for 
the phase diagram of doped manganites in terms of a 
model including the strong on-site repulsion, orbital de- 
generacy, and anisotropic covalencyEj. Based on this 
MFT, we first presents the spin wave dispersion in terms 
of the random-phase-approximation (RPA). This repro- 
duces qualitatively the x-dependence of the stiffness and 
the anisotropy due to the cross-over from superexchange 
interaction (SE) to DE. Especially for the doped region, 
only when the orbital configuration becomes d x i_ y i, i.e., 
x > 0.2, the DE becomes appreciable and the in-plane 
spin stiffness grows-jfaoidJy. Observed values of the in- 
plane spin stiffnessE3BEZl 
estimated value with cL.2 
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agree quantitatively with the 
- y i -orbital ordering. T, 
derstood in terms of the orbital liquid picture 
implies that the Jahn- Teller (JT) phonon has less impor- 
tance on the spin dynamics. 
We start with the Hamiltonian 

H = E *7/4 7 <w 
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where 7 [= a(d x 2_ y 2), b(d 3z 2 _ r 2)] specifies the orbital and 
the other notations are standard. The transfer integral 
tj- depends on the pair of orbitals (7,7') and the di- 
rection of the bond (i,j)§3. The spin operator for the 
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e g electron is defined as S, 
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the Pauli matrices a, while the orbital isospin operator 
is defined as Tj = | d^^S^td^i a . Jh is the Hund's 
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coupling between e g and t2 g spins, and Js is the AF 
coupling between nearest neighboring t2 g spins. H on s i te 
represents the on-site Coulomb interactions between e g 
electrons. Coulomb interactions induce both the spin 
and orbital isas-pjn moments, and actually H on s n e can 
be written asE 



E 



(pr. f + aS 2 egl ) 
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where the coefficients of thi 
i.e., a and /3, are given by^ 1 



and isospin operators, 
a = U - ir > , and 
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(3 = U — > . The parameters a, fi, to, used in the 
numerical calculation are chosen as to ~ 0.72 eV, U — 
6.3 eV, and—/ = 1.0 eV, being relevant to the actual 
manganites.Ej 

In the path-integral quantization, we introduce the 
Stratonovich-Hubbard fields s and 0t, representing the 
spin and orbital fluctuations, respectively. With the large 
values of the electron-electron interactions above, both 
(ps and (pT are almost fully polarized.EJ The MFT corre- 
sponds to the saddle point configuration of ips and (fix- 
We consider four kinds of spin alignment in the cubic 
cell: F-, A-, C- and G-type. As for the orbital degrees of 
freedom, we consider two sublattices /, and II, on p.ch 
of which the orbital is specified by the angle 6i t u asc3 
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We consider four types, i.e., F-, A-, C-, G-type also for 
the orbital ordering. Henceforth, we use a notation such 
as spin A. orbital G (^Jj^) etc - -*- n MFT, the most 



stable ordering is given by 

x = 0.0 Spin A Orbital C 

x = 0.1 Spin F Orbital C 

x = 0.2 - 0.4 Spin A Orbital F 

x = 0.5 - 0.9 Spin C Orbital F 



(60,- 
(80,- 
(0,0) 
(180,180). 
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As for x = 0, we further introduced the JT effectl 
putting the observed distortion of the MnOg octahedra.ES 
RPA corresponds to the Gaussian fluctuation around 
MFT, and the contribution to the spin wave effective 
action from the e 9 -electrons Ssw is obtained as the ex- 
pansion around the saddle point. 

<5sw = /] K * (?, SI) 7T {q s + q, SI) -it (-qs~q, -SI) 

+J2 K x(& n) *(qs + q, Sl)-{nXTr(-q s -q,-Sl)}. (4) 

where qs(= —qs) is the wavevector and ft (|n| = 1) is 
the direction of the ordered magnetic moment, and if 
is the fluctuation perpendicular to it. Because the spin 
wave is the Goldstone boson, the condition K n (0, 0) = 
, K x (0, 0) = 0, can be derived. Coefficient of the diago- 
nalized quadratic form is obtained as Ki(i) = K^HKx , 
zero-point of which (K^^ (q, SI = — iu>) = 0) gives the 
dispersion relation of the excitation lo = w(q). How- 
ever in this paper we focus on the static spin stiffness 
rather than the dynamic spin wave velocity because (a) 
at x — the spin stiffness is correctly reproduced to be 
of the order of J in the RPA while the spin wave velocity 
scales with t, and (b) for the metallic region, x ^ 0, the 
Landau-damping is not properly treated in our calcula- 
tion where the Brillouin zone is discretized and thus the 
gapless individual-excitation is not correctly evaluated. 
The static spin stiffness C a corresponds to the static re- 
sponse function for small |<f| as 
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and roughly reflects the exchange-interaction depending 
on x, where a = 1 (—1) corresponds to spin up (down), 
respectively. 

In RPA calculation the SE corresponds to the con- 
tribution from the inter-band transitions, while the DE 
from the intra-band ones. In this way, the present calcu- 
lation describes both SE and DE interactions, and hence 
their crossover in a unified way. Also the contribution 
from Js should be considered, the value of which is de- 
termined in the following way. We require that the exper- 
imentally observed anisotropy ratio of the spin stiffness 



R = y £' v j is reproduced when the calculated contri- 
butions from e s -electrons and that from are added. 
The observed value R = 7.6 for LaMnOaEj leads to an 
estimation as Js — 0.997 meV. As for AFM at x = 
the observed value R = 10.4 for Ndo.45Sr .55Mn03i 
gives Js = 1.4 mcV. These estimations are consistent 
rather with Js ~ 0.8 pfV estimated from the Neel tem- 
perature of CaMnOaEj than, the earlier meanfield es- 
timations Js ~ 8 meV.co Using these estimations, 
Js ~ 1 meV, we can estimate the spin wave stiffness for 
x = as ^total^total = -"-^ nicV, including the contri- 
bution from t 
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orbital. The corresponding expeumen- 
= 3.91 mcV in LaMn0 3 El with 
the reported lattice constants and the magnitude of spin 
moment, Stotai = 3/2 + 1/2(1 — x). The discrepancy 
may be attributed to the complex lattice deformations 
such as the Mn-0 bond-length (JT-type distortion) and 
the Mn-O-Mn band-angle (orthorhombic distortion) ob- 
served at 0L3, which can also be an origin of the 
anisotropycjEa. 




A-type Insulator 
(spin A, 

orbital F(60,-60» 



3 Dimensional 
F-type Metal 



quasi- 1 Dimensional 
C-type Metal 
(spin C, 

(spin A, orbital F<0,0) orbital F(180,-180)) 



2 Dimensional 
A-type Metal 



fspin F, orbital F(B0,-B0)) 



FIG. 1. Doping-dependence of the spin stiffness. The or- 
bital and the spin structure is optimized at each point. The 
enhancement of the spin-stiffness and the cross-over of the 
dimensionality are seen with increasing x. 
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We now turn to the doped case x ^ 0. Fig. [j] shows 
the doping-dependence of the total stiffness calculated for 
the optimized spin/orbital structure at each x. Firstly 
the spin stiffness due to the double exchange interaction 
scales roughly with x because the orbital is almost fully 
polarized, while in the absence of the orbital polarization 
it scales with the electron density (f — x) rather than the 
hole x for small x. The observed stiffness enhancement 
with increasing x even in the metallic regior£3 therefore 
also supports the large orbital polarization due to the 
strong Coulomb interactions. As x increases, the spin 
structure changes from spin A-type insulator at x = 
into the nearly isotropic FM, to the AFM with two- 
dimensional d x 2_ y 2 orbital alignment, and to the spin C 
metal with d 3z 2_ r 2 orbitalE3 Accordingly, the in-plane 
stiffness shows an increase, moderately at the beginning 
and then rapidly in the region of AFM. This reflects the 
fact that the DE is the most effective and prefers the 
d a; 2_j / 2-orbital, i.e., the DE is basically two-dimensional 
with the e g -orbitals. In the spin-C-metal for x > 0.4, 
one-dimensional orbital along (OOl)-direction gives rise 
to a steep increase of the stiffness in this direction. 

The observed anisotropy of the spin stiffness is deter- 
mined by the long range ordering of the orbitals. Fig. 
|l| also represents the cross-over of the dimensionality 
which— we proposed in the previous report .E2I Yoshizawa 
et alxl observed the reentrant of such two-dimensional 
anisotropy of the stiffness for Ndo.45Sro.55Mn03, be- 
ing consistent with our result. Quasi-one-dimensianal 
anisotropy is predicted for Ndi-^Sr^MnOs (x > 0.6)li3. 

The in-plane spin stiffness ^totai^totai m Fig- could 
be compared— with the experiments. In Lai-^Sr^MnOs, 
Endoh et alS3 observed the plateau of the velocity v x in 
the orbital-ordered insulating state up to x ~ 0.12 and 
then the velocity increases in the FM phase. Comparing 
this with the calculation above, it seems that the mod- 
erate increase up to x ~ 0.15 in Fig. [l] corresponds to 
the plateau, while the rapid increase fo«: > 0.15 to the 
increasing velocity observed by Endoh.ll3 Then orbital- 
ordered FM state in Fig. 1 corresponds to the insu- 
lating spin F phase in experiments. Both the FM and 
AFM phases in experiments, on the other hand, seems 
to corresponds to the AFM with d x 2_ y 2 orbital ordering 
in the calculation— This fits well orbital liquid picture 
by Ishihara et alx3, In a perfectly cubic system the or- 
bital state in FM is described as the resonance among 
d x 2_ y 2 1 d y 2_ z 2, and d z 2_ x 2. In the actual CMR-,com- 
pounds, however, the slight lattice distortiorJilta may 
breaks the cubic symmetry to stabilize d x 2_ y 2 though it 
is still accompanied with large fluctuation around it. 

Now we turn to the absolute value of the spin stiff- 
ness in FM phase. With the reported lattice constants 
the experimental values of the spin-stiffness, ^otai'^totai' 
are 11.61 meV fox La .7Sr . 3 MnO 3 y and 10.24 meV for 
Ndo.7Sr .3Mn03E3, respectively. These are in a good 
agreement with </ t otai^otai =10-53 meV estimated by 
RPA here with x = 0.3, d x 2 _ y 2 -orbital ordering (a sim- 



ple tight binding estimation with d 2 ,2_ J( 2-orbital also gives 
similar value). This agreement implies the large orbital 
polarization in FM phase with d x 2_ y 2 at least locally. 

This o rb i ta l-Liquid picture alsOj-explains the spin wave 
softeningc3ll3'E3 and spin cantingp3 observed in this sys- 
tem. Some theoretical works shows that the orbital fluc- 
tuation such an orbital liquid state leads to the softening 
of the spin wave dispersion-Jiear the zone boundarytll 
with the anisotropic featureEa (the. softening almost dis- 
appears along (-7T, 7r, 7r)-directional3'Ea). As for the spin 
canting, the observed canting in the metallic p^jegion 
(Ndo. 5 Sro. 5 Mn03) with the FM/AFM transition^ can- 
not be explained unless the planer orbital d x 2_ y 2 realizes 
in FM phasetil (observed slight lattice anisotropy can not 
stabilize such a planer orbital without the occurrence of 
the orbital liquid state). 

An important conclusion from the agreement between 
the experiments and RPA calculation of the stiffness con- 
stant is [that the polaron effect is small in the metal- 
lic stateEj at least on the spin dynamics. Polaron 
should reduce the DE interaction in the doped region 
via a bandwidth reduction by a factor of < xjXj >= 
cx P [-£ 9 Kl72] (u q = GtoHXe*** -e*«-*0), where 

X, t = exp[Y /q e lq ' Rl (g q /uj q )(b q - bL q )} is a factor en- 
countered in the canonical transformation eliminating 
the coupling between electrons and polaronic bosons, 
J2i,crJ2 q 9q(b q + bL^d^dur, with the coupling constant 
g q and phonon frequency uj q S3 On the other hand, for 
x = 0, the SE under the coupling with the polaron is 
given by, 

J = 4|^| 2 ^ /3 drG2(r)(xt(r)A,(r)Aj(0)X l (0)) , (6) 

where Gq(t) = e~ Ur ^ 2 is the Green's function for local- 
ized electrons. Because we are interested in the large 
U case, the integral is determined by the small t re- 
gion where ^(t^t)*] (0)A 2 (0)> S e" Ar (A = 
£ ? w g |u 9 | 2 ). Then the polaronic effect is to replace U 

by U + A in the expression for </, which is a minor cor- 
rection when U >> Aa, being in sharp contrast to DE 
discussed above. Polaronic effect should therefore correct 
the RPA-estimation of the stiffness-enhancement as x in- 
creases to be smaller. Agreement between the observed 
and estimated stiffness for DE implies therefore that the 
spin dynamics is not so affected by, the polaron. This is 
also pointed out by Quijada et alL3 

In summary, we have studied the role of orbitals in the 
spin dynamics of Ri^ x A x yhi0 3 . Comparing the experi- 
ments with the RPA calculation based on the mean field 
theory, we conclude the followings. (a) ^-dependence of 
the stiffness-enhancement suggests the large orbital po- 
larization, (b) the double-exchange interaction prefers 
d x 2_ y 2 orbital and is basically two-dimensional interac- 
tion, which leads to the large anisotropy of the spin dy- 
namics, (c) the agreement between experiments and RPA 
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results strongly suggests that the spin dynamics is not so 
affected by the JT polaron. 
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